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The process e + e~ — > 3 jets offers the opportunity to measure the strong coupling 
constant. For an accurate determination, precise theoretical calculations are nec- 
essary. I will give an overview on the status of the next-to-next-to-leading order 
calculations. 



1 Introduction 

Today's and future precision measurements in particle physics require accurate 
theoretical predictions in order to extract the values of the fundamental pa- 
rameters of the theory from experiment. The precise values of these constants 
also serve to constrain the parameter space of models of new physics. To this 
aim existing next-to-leading order (NLO) predictions have to be extended to 
next-to-next-to- leading order (NNLO). While for inclusive observables like the 
total hadronic cross section in e + e~ annihilation this step has been taken long 
time ago, the situation is much more complicated for exclusive quantities like 
event shapes or jet rates. One prominent process where a NNLO calculation for 
fully exclusive observables is desirable is e + e~ — > 3 jets. The strong coupling 
constant a s can be measured by using the data for e + e~ — > 3 jets. A NNLO 
calculation is expected to reduce significantly the theoretical uncertainty in 
the extraction of a s . A program for a fully exclusive NNLO calculation is 
flexible and allows to take into account complicated detector geometries and 
jet definitions. The observable may even be defined after the program has 
been written. The only requirement on the observable is infrared-safety. To 
construct an NNLO program the following ingredients are needed: 

- the relevant scattering amplitudes up to two loops; 

- a method to cancel infrared divergences; 

There has been significant progress in the past years for each of these items 
and I will review them in the remainder of this talk. 



2 The Calculation Of Two-Loop Amplitudes 

The NNLO calculation of e + e _ — > 3 jets requires the following amplitudes: the 
Born amplitudes for e + e~ — > 5 partons, the one-loop amplitudes for e + e~ — > 
4 partons and the two-loop amplitudes for e + e _ — > 3 partons. The necessary 
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Born amplitudes 1,2 and the one-loop amplitudes 3,4 have been obtained already 
a while ago. The most complicated parts are the two-loop amplitudes. This is 
due to new types of Feynman integrals, corresponding to Feynman diagrams 
with two loops and an external off-shell leg. Methods invented to tackle these 
integrals comprise the use of the Mellin-Barnes formula 5 , the application of 
differential equations and integration-by-parts identities 6 as well as the use of 
nested sums 7 . The results are expressed in terms of (multiple) polylogarithms, 
thus extending the well known set of basic functions for one-loop integrals 
(logarithm and dilogarithm) and making contact with recent developments 
in mathematics. In addition to the results for the basic integrals, various 
reduction algorithms are also used: An algorithm by Tarasov 8 allows to convert 
integrals with the loop momentum in the numerator into scalar integrals with 
possibly raised powers of the propagators and shifted dimension. An efficient 
algorithm by Laporta 9 expresses a large set of unknown integrals in terms 
of a few "master" integrals. With these tools the two-loop amplitudes for 
e+e~ — ► qgq have been calculated 10,11 . The two-loop amplitude has the colour 
decomposition 

4 2) = iv 2 4 2 l+4:i + ^4:i+^4!l + f 45+^ 2 45 

+A {2) + A (2) 

(2) (2) 

A 3 vec and .A3 ax are the contributions resulting from the vector and axial 

vector coupling of the photon/Z-boson to a closed quark loop. The Durham 

(2) 12) (2) 
group calculated the partial amplitudes A 3 [ - A 3 6 as well as A 3 vec . Our group 

(2) (2) (2) 

obtained results for A 3 4 - A 3 q. The results from the two groups for A 3 4 - 

(2) (2) 
A 3 6 agree analytically. The partial amplitude A 3 ax remains to be calculated, 

although it is expected that this partial amplitude gives a negligible numerical 

contribution. 

3 Cancellation Of Infrared Divergences 

At the next-to-next-to-leading order level the ingredients for the third order 
term in the perturbative expansion for quantities depending on n resolved 
"hard" partons are the already mentioned n-parton two-loop amplitudes, the 
(n + l)-parton one-loop amplitudes and the {n + 2) Born amplitudes. Taken 
separately, each one of these contributions is infrared divergent. Only the sum 
of all contributions is infrared finite. Infrared divergences occur already at 
next-to-leading order. At NLO real and virtual corrections contribute. The 
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virtual corrections contain the loop integrals and can have, in addition to ul- 
traviolet divergences, infrared divergences. If loop amplitudes are calculated in 
dimensional regularisation, the IR divergences manifest themselves as explicit 
poles in the dimensional regularisation parameter e = 2 — D/2. These poles 
cancel with similar poles arising from amplitudes with additional partons but 
less internal loops, when integrated over phase space regions where two (or 
more) partons become "close" to each other. In general, the Kinoshita-Lee- 
Nauenberg theorem 12,13 guarantees that any infrared-safe observable, when 
summed over all states degenerate according to some resolution criteria, will 
be finite. However, the cancellation occurs only after the integration over the 
unresolved phase space has been performed and prevents thus a naive Monte 
Carlo approach for a fully exclusive calculation. It is therefore necessary to 
cancel first analytically all infrared divergences and to use Monte Carlo meth- 
ods only after this step has been performed. At NLO, general methods to 
circumvent this problem are known. This is possible due to the universality of 
the singular behaviour of the amplitudes in soft and collinear limits. Exam- 
ples are the phase-space slicing method and the subtraction method. Within 
the subtraction method 14,15 one subtracts a suitable approximation term da A 
from the real corrections da R . This approximation term must have the same 
singularity structure as the real corrections. If in addition the approximation 
term is simple enough, such that it can be integrated analytically over a one- 
parton subspace, then the result can be added back to the virtual corrections 
d<r v . Since by definition da A has the same singular behaviour as d<r R , the 
combination (da R — da A ) is integrable and can be evaluated numerically. Sec- 
ondly, the analytic integration of da A over the one-parton subspace will yield 
the explicit poles in e needed to cancel the corresponding poles in da v . At 
NNLO this generalizes as follows 16 : 

<0>r L ° = / O n+2 da n % O n+1 o - O n o da^ 

+ J O n+1 d<r n % + O n+1 o da£# - O n o rial 1,1 ) 
+ J On d*W + O n o da<£& + On ° da™. 

Here da^+i is a subtraction term for single unresolved configurations of Born 
amplitudes. This term is already known from NLO calculations. The term 
da>n' 2 ^ is a subtraction term for double unresolved configurations. Finally, 
dan' 1 ^ is a subtraction term for single unresolved configurations involving one- 
loop amplitudes. 
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To construct these terms the universal factorisation properties of QCD 
amplitudes in unresolved limits are essential. QCD amplitudes factorise if they 
are decomposed into primitive amplitudes. Primitive amplitudes are defined 
by a fixed cyclic ordering of the QCD partons, a definite routing of the external 
fermion lines through the diagram and the particle content circulating in the 
loop. One-loop amplitudes factorise in single unresolved limits as 17 > 18 ' 19 > 20 

4 1} = Sing^.^ + SW 1 - 1 ).^ ),. 
Tree amplitudes factorise in the double unresolved limits as 21 < 22 ,23,24,25 

4° } = Sing^).4 2 2 . 

The subtraction terms can be derived by working in the axial gauge. In this 
gauge only diagrams where the emission occurs from external lines are relevant 
for the subtraction terms. Alternatively, they can be obtained from off-shell 
currents and antenna factorisation 26 . 

4 Outlook 

In this talk I reviewed the status of NNLO 3-jet calculations. With the progress 
we witnessed in the field in the last years we can expect to obtain numerical 
results rather soon and to extend existing numerical programs for NLO predic- 
tions on e+e~ — ► 4 jets 27 - 28 - 29 ' 30 towards NNLO predictions for e+e~ — > 3 jets. 
In fact, a result for a particular colour structure for the thrust observable was 
announced recently 31 . 
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